A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 

02.02.1; 11.09.3; 11.10.1; 11.03.4: Coma;12.12.1 



ASTRONOMY 

AND 
ASTROPHYSICS 



66 Black Hole Energy Release to the Gaseous Universe 

o\ ■ 

^ ■ Torsten A. Enfilin 1 , Yiping Wang 12 ' 3 , Biman B. Nath 4 , and Peter L. Bi 



siermann 



a ■ 



Max-Planck-Institut far Radioastronomie, Auf dem Hiigel 69, D-53121 Bonn, Germany 
Purple Mountain Observatory, Academia Sinica, Nanjing 210008, China 
Bergische Universitat Wuppertal, D-42097 Wuppertal, Germany 
Raman Research Institute, Bangalore 560080, India 



Received ??? , Accepted 



??? 



> 

in : 
o , 

m ■ 

o : 

00 , 
On ■ 

^ : 

Oh; 

6 ' 
a: 



Abstract. We estimate the energy release of black hole 
formation to the intra-cluster medium of the Coma clus- 
ter and find that this is comparable to the present day 
energy content. Therefore an energetic and maybe hydro- 
static influence is possible. Our calculations rely on the 
assumption of an universal black hole to galaxy mass ratio 
(more exactly: spheroidal mass component of the stellar 
population), for which there is growing evidence. On a cos- 
mological scale, there is also an energy release of black hole 
formation comparable to what is expected to be present 
within the thermal gas, caused by the process of structure 
formation. This indicates an important dynamical influ- 
ence, neglected by present day structure formation sim- 
ulations. This estimate of cosmological black hole energy 
release is independent of the black hole to galaxy mass 
ratio, but consistent with its value. 
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host galaxies energetically and dynamically. In this Let- 
ter we present a comparison of the BH energy release and 
the thermal energy of the quasar environment in order to 
demonstrate the possible importance of this influence. The 
efficiency of energy dissipation of the accreting matter to 
a Schwarzschild BH is e w 0.06, whereas to a maximal ro- 
tating BH e w 0.3 (Thome 1974; Laor & Netzer 1989). We 
assume £ii gri t ~ 0.1 to be the mass-to-light conversion rate, 
and £th+nth ~ 0.1 for the production efficiency of thermal 
and nonthermal energy release in relativistic particles and 
magnetic fields. Therefore we use an intermediate value of 
e w 0.2. We adopt a conservative value of ?/bh ~ 0.002, 
since we like to include the fraction of galaxies containing 
massive BHs into this number. If rybh ~ 0.006 our esti- 
mate of energy release to the Coma cluster is still correct 
if only one third of all galaxies contain massive BHs. But 
the cosmological estimate is not affected, since it depends 
only on the ratio eth+nth/^iight and the observed quasar 
light. Further, we use H = 50/i5okms _1 Mpc -1 , f2 Q = 1, 
A = 0, and indicate the scaling of all quantities with the 
Hubble constant. 



1. Introduction 

There is growing evidence for the existence of black holes 
(BHs) in all galaxies, being remnants of the strong quasar 
activity in the young Universe (Lynden-Bell 1969; Rees 
1989; Haehnelt & Rees 1993; van der Marel 1997; Ford 
et al. 1997; Hasinger 1998; and references therein). Esti- 
mates of the ratio of the BH mass to that of the sphe- 
roidal component of the stellar population of the host 
galaxy seem to converge to r/hh ~ 0.002 — 0.006 (Kor- 
mendy & Richstone 1995; Faber et al. 1997; Magorrian et 
al. 1997; Silk & Rees 1998; Wang & Biermann 1998). This 
allows to estimate the total BH mass within a popula- 
tion of galaxies. The former BH growth was accompanied 
by an output of large amounts of energy by radiation, 
heating of the ambient medium and outflows of relativis- 



2. The Coma Cluster 

First, we estimate the thermal energy content of the 
Coma cluster. The electron density of the intra-cluster 
medium of the Coma cluster can be described by a [3- 
model: n c (r) = n . o [l + (r/r corc ) 2 ]~ 3/3t =/ 2 ; where n Cj0 = 

3 ■ 10" 3 

^50 2 Cm 3 1S tne central electron density, r coro — 
400/i 50 1 kpc the core radius, and (3 C = 0.75 gives the 
slope (Briel et al. 1992). The gas extends to a radius of 
R cluster ~ 5 h 5 Q Mpc, where an accretion shock of the 
infalling matter marks the cluster boundary (EnBlin et 
al. 1998). The integrated gas mass within this radius is 
4.2 • 10 14 h 5 Q^ 2 Mq. The central temperature of the clus- 
ter is 8.2 keV (Briel et al. 1992), but there is evidence 
for a temperature decrease from the center to the bound- 
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with (3kT = 0.09, leading to a temperature drop by two 
within 5 Mpc. The resulting thermal energy content 
within -Rciustor is 1.3 • 10 64 h^ /2 erg. 

This has to be compared with the energy release of 
BH formation. Since the galaxies in Coma are mainly el- 
lipticals, only their total mass has to be estimated and 
multiplied with r/bh hi order to get the mass of BHs in 
Coma. The mass of the ellipticals can be derived from 
their luminosity function and their mass to light ratio. 
The R-band luminosity function of early type galaxies in 
the central 700 arcmin 2 of Coma is given in Seeker & Har- 
ris (1996). The complete luminosity function can be es- 
timated by assuming these galaxies to be within a cylin- 
drical volume, and correcting to the whole cluster vol- 
ume, with the help of the spatial distribution of galax- 
ies. The radial distribution of galaxies within Coma is 
given by n ga i ~ (1 + (V/r gal ) 2 )~ Q sai ; w ith a ga i = 0.8 and 
r sa i = 160/igo kpc (Girardi et al. 1995). Including into 
this the luminosities of NGC 4874 and NGC 4889 (Strom 
& Strom 1978), which were not included into this lumi- 
nosity function, and applying a correction of V-R = 0.85, 
which is typical for ellipticals (Tinsley & Gunn 1976), we 
get the V-band luminosity function. The galactic mass 
can now be integrated with the help of the mass to (V- 
band) light ratio given by Magorrian et al. (1997). It 
is M gal = 2.9 • 10 13 ft,^ M , giving a total BH mass of 
M bh = 5.8 • 10 10 (?7 bh /0.002) M . The gravitational 
energy, which was freed during the BH formation, is there- 
fore E hh = sM hh c 2 = 2.1 • 10 64 h^ 1 (e/0.2) (?7 bh /0.002) erg. 
The thermal plus nonthermal energy release (one half of 
the total) is therefore comparable to the present day ther- 
mal energy content of Coma within a 5 Mpc radius. 

A different approach to the nonthermal energy release 
of BH formation to clusters of galaxies was taken by Enfilin 
et al. (1997) by an estimate of the jet-power of radio 
galaxies, integrated over cosmological epochs. A compar- 
ison with the typical thermal energy content of cluster 
showed, that both numbers could be comparable on clus- 
ter scale, if a considerable amount of energy is ejected 
by the not radio-emitting and therefore invisible relativis- 
tic protons. Our above estimate of the nonthermal en- 
ergy release is completely independent of this calculation, 
and therefore demands such energetic protons, having en- 
ergy densities much higher than the electrons within the 
relativistic plasma flowing out of radio galaxies, consis- 
tent with theories of the origin of the observed ultra-high- 
energy cosmic rays (Rachen & Bicrmann 1993; Bicrmann 
1997) and recently discovered TeV 7-rays from blazars 
(Mannheim 1998). 

Since also the nonthermal energy should be stored in 
the intra-cluster medium for at least a Hubble time, we 
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cluster), Enfilin et al. (1997), and Berezinsky et al. (1997). 
We note, that details of the injection of relativistic plasma 
from quasars into the intra-cluster medium and also the 
later evolution of this plasma need further investigations. 
E.g. in speculative scenarios of dwarf galaxy formation in- 
duced by violent expansion of quasar outflows in the early 
Universe large fractions of the released energy are trans- 
formed to heat in the ambient medium by shocks and adi- 
abatic compression (Natarajan et al. 1997; Silk & Rees 
1998; and references therein). The question if the non- 
thermal plasma influences the hydrostatics of the cluster, 
and should be taken into account into (uncertainties of) 
X-ray mass determinations, can not be answered yet, since 
this depends on the amount and also crucial on the spatial 
distribution of the nonthermal phases. 

At the location of the radio outflow from Perseus A 
an X-ray hole is present in the thermal emission of the 
Perseus cluster (Bohringer et al. 1993), demonstrating 
that the radio plasma rather displaces the thermal gas, 
driven in this case with a power of probably more than 
10 46 ergs~ 1 (Heinz et al. 1998), than mixes with it. Thus 
the radio plasma might remain within the intra-cluster 
medium, invisible since the electrons cool down in a cos- 
mologically short time, and might be reactivated when 
a cluster merger event reaccelerates the electron popu- 
lation, as it could be the case for the 3 h§ Mpc sized, 
diffuse radio halo of Coma (Deiss et al. 1997; and refer- 
ences therein). Also the peripherally located cluster ra- 
dio relics (e.g. 1253+275 in Coma) can be understood 
in terms of old reactivated radio plasma, having passed 
through a cluster accretion- or merger shock wave (Enfilin 
et al. 1998). Evidence for the existence for a necessary, 
old, but large relativistic electron population within the 
intra-cluster medium, having energies below what is visi- 
ble in the radio, can be seen in the recently detected EUV 
excess of the Coma cluster (Lieu et al. 1996), if this is 
explained by inverse-Compton scattered microwave back- 
ground photons (Hwang 1997; Enfilin & Biermann 1998; 
Sarazin & Lieu 1998). 



3. The Universe 

It is also possible to compare the BH formation energy 
release to that of the thermal gas on a cosmological scale. 
The integrated energy density in quasar light emitted 
during (and estimated for) the epoch of quasar activity 
is (0.9 - 1.3) • 10~ 15 ergcm- 3 (Chokshi & Turner 1992; 
an earlier, more conservative estimate can be found in 
Soltan 1982). This number is given per comoving vol- 
ume, and it is independent of cosmology. Chokshi & 
Turner (1992) assumed an efficiency of light production 
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A large fraction of the energy output is injected into 
clusters, roughly / d = 0.3 - 0.5 (Enfilin ct al. 1997), 
which contain less than 10% of all baryons, and there- 
fore get more energy per baryon from BHs than the gas 
outside clusters. Distributing 60% of the BH energy re- 
lease onto the 90% of the gas which is outside clusters 
gives an average temperature of fcT(z; n j) < (0.9 — 1.3) keV 
(eth+„th/eiight) ((1 - /d)/0.6) (ft b /i§ /0.05 J" 1 at the red- 
shift of injection. This estimate assumes complete dissi- 
pation of the nonthermal energy in magnetic fields and 
relativistic particles. This is not realistic, but gives us up- 
per limits to the heating and Comptonization parameter. 
The true thermal heating might only be some fraction of 
this number. 

In a homogeneous universe comoving energy densities 
have to be redshifted by (1 + Zi n j)~" due to adiabatic ex- 
pansion losses, with a — 3(7 — 1) and 7 the adiabatic 
index. This correction has to be applied for the quasar 
light, which has 7ii g ht = 4/3 and therefore an g ht = 1- The 
average injection redshift of the luminosity function used 
in Chokshi & Turner (1992) is z; n j = 2.34, giving a present 
day quasar light density of (0.3 — 0.4) • 10~ 15 erg cm -3 . The 
redshift correction for BH thermal energy release should 
be different, due to 7 gas = 5/3, but also since this heat is 
injected into the environment of the quasar host galaxies, 
which is the inter-galactic medium of the filaments and 
sheets of the large scale structure and therefore partly de- 
coupled from the Hubble flow. We assume that the den- 
sity in sheets and filaments decreases with (1 + z) d only 
due to the action of the Hubble flow in the direction of 
these structures of dimension d: we adopt d = 2 for sheets 
and d = 1 for filaments. Since most galaxies are located 
in filaments, we use d = 1, but give the scaling for dif- 
ferent values. This is conservative, since these structures 
are forming, and the ongoing infall of matter could lead 
to compression, whereas the Universe still expands. Thus, 
a% as = d(j gas — 1) = 2/3 gives a present day average heat 
(from BHs) outside clusters of 

kT(z = 0) < (0.4 - 0.6) keV --±-- -— ^~ 

£light 0.6 

3.34- 2 / 3 V 0.05 J ' 

(independent of cosmology) again assuming complete dis- 
sipation of nonthermal phases (therefore <). The corre- 
sponding Comptonization parameter y < (0.7 — 1.2) • 
lO~ 5 h§Q, estimated for this adiabatic cooling history and 
an Einstein-de-Sitter Universe (using Q Q = increases 
y by only 50%), does not violate the present day up- 
per limit of y < 1.5 • 10" 5 (Fixscn ct al. 1996). Addi- 
tional to the dependence on /150, y scales with (1 + Zi n j) 3 / 2 
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The true cosmological BH heating of gas in cosmolog- 
ical filaments might be lower than this heat, due to the 
nonnegligible fraction of nonthermal energy which does 
not dissipate, and which is mainly invisible to Comp- 
tonization measurements. But this phase cools slower than 
the thermal gas, due to a smaller adiabatic index. The 
present day BH energy release should therefore be slightly 
higher than the thermal energy given above (present day 
BH released thermal and nonthermal energy outside clus- 
ters per cosmological volume: (2.5 — 4.0) • 10 -16 erg cm -3 ). 
This is large enough in order to be energetically and dy- 
namically important. For comparison: numerical simula- 
tion for the evolution for IGM assumes only photoioniza- 
tion from quasars and get a temperature increase of sev- 
eral eV (e.g. Miralda-Escude et al. 1996; and references 
therein). Simulations of structure formation predict tem- 
peratures of 0.1 — 0.5 keV within filaments, resulting from 
shock heating at the accretion shocks of filaments (Kang 
et al. 1996). This demonstrates the need for simulations 
of structure formation, which take into account the back 
reaction of galaxy formation induced BH growth onto the 
inter-galactic medium. 

4. Conclusions 

We estimate the contribution of energy release of black 
hole formation during the quasar active epoch of the Uni- 
verse for the energy budget of the Coma cluster of galax- 
ies and also on a cosmological scales in order to inves- 
tigate their possible energetical and dynamical influence. 
Although uncertainties are large, we find, that these en- 
ergies are comparable to that of the thermal gas of the 
environment. In the case of galaxy clusters a possible hy- 
drostatic influence of nonthermal phases was proposed as a 
possible systematic effect for mass determinations of clus- 
ters, and therefore for f2 (Loeb & Mao 1994; Enfilin et al. 
1997). This Letter shows that sufficient energy seems to be 
injected (depending on the black hole efficiency and black 
hole to spheroidal mass ratio), but the unknown details 
of thermalization and spatial distribution of nonthermal 
phases within the intra-cluster medium do not allow a final 
answer to this question. Since large amounts of magnetized 
plasma should have flown out of quasars, this scenario also 
can help to explain the origin of the observed intra-cluster 
magnetic fields. On larger scales, the black hole energy re- 
lease can influence the large-scale structure formation by 
the feed back reaction during galaxy formation and black 
hole growth. This cosmological estimate is quite robust 
against changes in the unknown parameters of black hole 
energy release, since it depends only on the ratio of the 
thermal plus nonthermal to light production efficiency. 
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